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HIGHLIGHTS 


►  2D  models  with  randomly  oriented  microstructures  and  various  grain  sizes  are  employed. 

►  The  grain  boundary  (GB)  is  modeled  as  a  thin  layer  between  grains. 

►  The  relative  orientation  angle  between  neighboring  grains  affects  Li  diffusivity. 

►  The  diffusivity  of  a  small-grained  structure  is  determined  by  the  properties  of  the  GB. 

►  The  influences  of  the  GB  properties  and  the  orientation  angle  can  be  evaluated  separately. 
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Layered  Li  transition  metal  oxides  are  widely  used  as  active  materials  for  the  positive  electrode  of  Li-ion 
rechargeable  batteries,  where  intercalation  of  Li  in  the  metal  oxide  is  a  fundamental  phenomenon  that 
determines  the  performance  of  the  batteries.  The  intercalation  process  is  significantly  affected  by  the 
crystal  anisotropy  and  grain  boundaries,  particularly  for  all-solid-state  thin  film  batteries.  Therefore,  to 
improve  the  batteries,  a  thorough  understanding  of  the  intercalation  process  on  the  nanometer  length 
scale  is  essential.  To  this  end,  we  have  proposed  phase-field  models  for  calculating  the  relation  between 
the  realistic  polycrystalline  microstructure  and  the  apparent  Li  diffusion  coefficient.  A  crystallographic 
orientation  was  randomly  allocated  to  each  crystal  grain  in  a  two-dimensional  polycrystalline  micro¬ 
structure.  The  simulation  results  show  that  the  apparent  Li  diffusivity  is  sensitive  to  the  diffusivity  of  the 
grain  boundaries,  the  spatial  distribution  of  the  crystal  orientation  for  each  grain,  and  the  grain  size.  The 
diffusivity  of  a  small-grained  structure  is  determined  by  the  properties  of  the  grain  boundary.  On  the 
other  hand,  the  diffusivity  of  a  large-grained  structure  depends  considerably  on  the  relative  orientation 
angle  between  neighboring  grains,  even  when  the  Li  conductivity  of  the  grain  boundary  is  large. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Layered  Li  transition  metal  oxides  are  widely  used  as  active 
materials  for  Li  ion  rechargeable  batteries  [1].  The  extraction  of  Li 
ions  from  and  the  insertion  of  Li  ions  into  a  layered  cobalt  dioxide 
matrix  [2-4]  are  fundamental  phenomena  that  determine  the 
battery  performance  [5].  This  structure  provides  two-dimensional 
(2D)  Li  diffusion  paths,  which  allows  topotactic  electrochemical 
reactions  to  occur  for  high-power  applications.  The  self-diffusion 
coefficient  of  Li  in  LixCo02  (Dseif)  has  been  investigated  using  both 
theoretical  and  experimental  approaches  [6-8].  The  value  of  Dseif 
estimated  from  muon-spin  spectroscopy  [6]  is  in  the  range  from 
1  x  10-10  to  1  x  10-9  cm2  s-1  for  LixCo02  (0.5  <  x  <  0.8)  at 


*  Corresponding  author.  Toyota  Central  R&D  Labs.,  Inc.,  41-1  Yokomichi, 
Nagakute,  Aichi  480-1192,  Japan.  Tel.:  +81  561  71  7127;  fax:  +81  561  63  5258. 
E-mail  address:  el044@mosk.tytlabs.co.jp  (S.  Yamakawa). 

0378-7753 /$  -  see  front  matter  ©  2012  Elsevier  B.V.  All  rights  reserved. 
http://dx.doi.Org/10.1016/j.jpowsour.2012.09.055 


a  temperature  of  300  K,  which  is  in  good  agreement  with  the  value 
calculated  using  a  first-principles  method  [8].  On  the  other  hand, 
the  chemical  diffusion  coefficient  (DChem),  which  is  evaluated 
by  electrochemical  measurements,  varies  from  1  x  1CT12  to 
1  x  10-1°  cm2  s-1  [9,10]  for  (104)-textured  thin  films,  so  that  grains 
with  their  (104)  planes  parallel  to  the  Li  transport  direction  facili¬ 
tate  the  Li  transport  [11  ].  It  is  known  that  DChem  can  be  expressed  as 
the  product  of  Dseif  and  a  thermodynamic  factor  (w).  Since  the 
binary  solution  of  Li  vacancies  in  the  2D  layer  is  not  ideal,  the  value 
of  to  must  be  larger  than  1  for  a  wide  range  of  Li  concentrations. 
Therefore,  DChem  is  expected  to  be  larger  than  Dseif,  which  contra¬ 
dicts  the  experimental  data.  The  transport  properties  of  Li  seem  to 
be  limited  by  microstructural  inhomogeneities  such  as  the  crystal 
anisotropy,  grain  size,  and  grain  boundaries.  In  a  commercial  Li-ion 
battery,  the  electrolyte  that  intrudes  into  pores  within  a  secondary 
particle  may  provide  short-circuit  transport  paths  for  Li.  Never¬ 
theless,  experimental  observation  [12]  has  shown  that  secondary 
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particle  contains  numerous  grain  boundaries,  which  means  that 
a  primary  particle  within  a  secondary  particle  is  not  fully  isolated 
and  the  surface  does  not  fully  expose  to  the  electrolyte.  On  the 
other  hand,  in  the  case  of  the  all-solid-state  battery  with  a  flat  thin- 
film  electrode,  since  most  crystal  grains  do  not  directly  touch  solid 
electrolyte  particles,  the  inter-granular  diffusion  between  crystal 
grains  plays  a  crucial  role  in  the  Li  transport.  Therefore,  to  improve 
the  batteries,  a  thorough  understanding  of  the  intercalation  process 
on  the  nanometer  length  scale  is  essential.  This  paper  focuses  on 
gaining  a  quantitative  understanding  of  the  relation  between  the 
morphology  of  the  microstructure  and  the  Li  diffusivity  using 
numerical  simulation  techniques. 

Most  previously  reported  numerical  simulations  [13-15]  have 
adopted  several  assumptions  about  the  anisotropic  Li  intercalation 
process,  such  as  an  isotropic  grain  model  and  linear  Fickian  diffu¬ 
sion  of  Li  inside  the  grains.  In  these  models,  because  the  influence  of 
the  polycrystalline  anisotropic  microstructure  on  the  Li  mobility 
was  taken  into  account  through  apparent  diffusion  coefficients,  the 
reaction  and  the  diffusion  process  of  the  practical  battery  electrode 
were  greatly  simplified.  Consequently,  such  models  could  fail  to 
capture  important  aspects  of  the  charging  or  discharging  dynamics 
in  an  essentially  anisotropic  in  homogeneous  medium.  It  can  easily 
be  imagined  that  a  linked  approach  for  simultaneously  modeling 
the  microstructure  and  the  transport  will  be  indispensable,  because 
microstructural  modeling  that  does  not  consider  transport  cannot 
identify  optimal  transport  properties  for  a  particular  application. 
The  phase-field  method  has  recently  attracted  increasing  attention 
as  a  promising  technique  [16,17]  for  calculating  the  morphological 
characteristics  through  microstructural  evolution  on  a  nanometer 
length  scale.  The  method  has  also  been  applied  to  describe  Li 
diffusion  in  a  phase-separated  system  [18,19]  and  to  simulate  the 
anisotropic  electrochemical  strain  microscopy  (ESM)  response  of 
polycrystalline  LiCo02  [20].  In  this  paper,  an  integrated  computa¬ 
tional  approach  for  simulating  the  morphology  of  the  poly¬ 
crystalline  microstructure  and  evaluating  the  corresponding 
apparent  Li  diffusion  coefficient  (Dapp)  in  the  microstructure  is 
proposed.  A  realistic  2D  polycrystalline  model  has  been  generated 
from  the  phase-field  simulations.  A  crystallographic  orientation 
was  randomly  allocated  to  each  crystal  grain.  In  this  model,  the 
grain  boundary  (GB)  was  treated  as  a  thin  layer  between  two  crystal 
grains.  Dapp  was  obtained  from  the  temporal  evolution  of  the  Li 
concentration  by  using  the  Cahn-Hilliard  diffusion  equation  [21], 
which  treats  non-Fickian  Li  diffusion  in  combination  with  the 
anisotropic  self-diffusion  tensor.  A  phenomenological  constitutive 
relation  between  Li  conductivity  and  quantitative  microstructural 
features  such  as  average  grain  size  and  crystallographic  orientation 
was  directly  provided  from  the  simulation  results. 

The  outline  of  this  paper  is  as  follows.  In  Sec.  2,  the  main  char¬ 
acteristics  of  the  simulation  methods  used  to  model  the  poly¬ 
crystalline  microstructure,  the  Li  diffusion  properties,  and  their 
integration  are  described.  In  Sec.  3,  practical  applications  are 
illustrated  using  a  2D  system.  The  dependence  of  Dapp  on  the 
various  grain  sizes  and  crystal  orientation  angles  between  neigh¬ 
boring  grains  in  the  microstructure  models  is  discussed.  We  also 
used  the  phase-field  model  to  study  the  Li  segregation  in  poly¬ 
crystalline  Li*Co02  thin  film  during  the  discharge  process  at 
a  constant  current.  The  conclusions  are  given  in  Sec.  4. 


2.  Calculation  method 

2.1.  Cahn-Hilliard  diffusion  equation 

In  order  to  represent  the  isotropic  Li  diffusion,  we  introduce  the 
Cahn— Hilliard  diffusion  equation  expressed  in  the  following  form: 


c(l  c)Dseif 
RT 


a) 


where  c  is  the  local  concentration  of  Li  normalized  by  the 
maximum  Li  concentration  (cJJ0*),  so  that  it  is  nondimensional  and 
only  takes  values  between  0  and  1.  Dseif  is  the  self-diffusion  coef¬ 
ficient  of  Li.  R  and  T  are  the  gas  constant  and  the  absolute 
temperature,  and  /jl  represents  the  diffusion  potential  of  Li,  which  is 
defined  as  the  functional  derivative  of  the  total  chemical  free 
energy  (Gsys)  with  respect  to  local  concentration: 


^Gsvs 


(2) 


Gsys  is  expressed  as  the  volume  integral  of  the  local  chemical 
free  energy  (Gm)  and  the  gradient  energy  term  over  the  entire 
system: 


Jsys 


V 


j  [Gm  +  |(Vc)2]  dV, 


(3) 


where  Gm  is  the  Gibbs  energy  of  the  uniform  binary  solution  in  the 
Li-vacancy  system  and  k  is  the  gradient  energy  coefficient,  which  is 
estimated  from  the  interaction  parameters  obtained  by  fitting  the 
Gibbs  energy  function  [22]  to  the  mixing  enthalpy  of  the  disordered 
solution.  We  employed  the  mixing  enthalpy  value  estimated  from 
the  first-principles  calculation  given  in  Refs.  [23,24].  By  substituting 
Eqs.  (2)  and  (3)  into  Eq.  (1),  it  can  be  rewritten  as 


c(l  -  c)Dself 

RT 


Eh_kV2 

Oc 


(4) 


here,  0Gm/0c  is  evaluated  by  utilizing  the  electrochemical  potential: 
^  =  -nLiF(ULiCo02  -  l/ave),  (5) 

where  ny  is  the  number  of  electrons  in  the  charge— transfer  reac¬ 
tion,  i.e.,  ny  =  1  in  Li+  +  e_  =  Li,  and  F  is  Faraday’s  constant.  DLiCo02  is 
the  equilibrium  potential  of  LiCoC^,  which  is  evaluated  by  using  Eq. 
(6)  in  Ref.  [25] 
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LiCo  02 


(6) 


It  is  noted  that  DLiCo02  is  given  as  a  function  of  c.  The  parameters 
Cmax.i,  ?,  and  U f  were  determined  by  fitting  the  function  to  exper¬ 
imental  open-circuit  potential  data  [25].  The  average  potential 
(Dave)  is  given  by  the  expression 


L^ave  =  (Guco02  “  GCo02  ~  GIi)/(_n Li*7)-  (7) 

To  verify  the  accuracy  of  the  present  phase-field  model,  we  used 
it  to  estimate  the  thermodynamic  factor  (w)  defined  as  the  ratio  of 
Dc hem  to  the  given  Dseif  in  homogeneous  LiCo02.  A  one-dimensional 
model  with  an  electrode  thickness  of  1  pm  was  employed  to  mimic 
the  potentiostatic  intermittent  titration  technique  (PITT)  [26].  The 
temperature  was  set  at  300.15  K.  A  potential  step  of  2  mV  was 
applied,  and  the  current  as  a  function  of  time  was  recorded.  DChem 
was  calculated  from  the  time  dependence  of  the  current.  The 
accuracy  of  this  simulation  program  is  verified  through  the 
comparison  of  the  calculated  oj  with  the  experimental  result  [10],  as 
shown  in  Fig.  1.  It  can  be  seen  that  w  exhibits  considerable  varia¬ 
tions  near  the  composition  of  <5  =  0.5  (Lio.sCoC^),  which  are  induced 
by  an  increase  in  the  degree  of  Li-vacancy  ordering  [27].  The 
simulated  result  was  overall  in  good  agreement  with  the  measured 
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Co02) 


where  the  coefficient  (3  is  introduced  to  represent  that  the  Li 
diffusivity  along  the  GB  core  and  that  across  the  interface  between 
a  grain  and  a  GB  layer  are  different  from  the  diffusivity  within  the 
grain: 


9c  _  _2_  f  j  j  \  |  d  C(1  Qffi^self 

at  ~  h  \Ji=h/2  Ji=~h/2)  +  as  rt  as 


(10) 


here,  h  is  the  GB  width  Ji=h/2  or  -h/2  denotes  the  flux  of  Li  across  the 
interface  between  a  grain  and  a  GB  layer,  and  s  is  the  direction  along 
the  GB.  Calculations  were  performed  for  such  an  in  homogeneous 
system  consisting  of  the  grain  and  the  GB  by  solving  the  two  types 
of  differential  equations,  Eqs.  (8)  and  (10),  simultaneously. 


2.3.  Polycrystalline  microstructure 


Fig.  1.  Variation  of  the  thermodynamic  factor  of  Li,  which  is  defined  as  the  ratio  of  the 
chemical  diffusion  coefficient  to  the  self-diffusion  coefficient,  with  composition. 

result  [10],  legitimating  the  modeling  and  the  parameters  used  for 
the  simulation. 

2.2.  Li  diffusion  in  poly  crystalline  material 


In  order  to  quantitatively  discuss  the  relation  between  the 
polycrystalline  microstructure  and  the  Li  diffusivity,  we  prepared 
various  types  of  polycrystalline  microstructures,  as  shown  in  Fig.  3. 
The  microstructures  were  obtained  by  using  the  multi-phase-field 
algorithm  [29]  under  periodic  boundary  conditions.  The  ideal 
grain  growth  for  a  system  was  simulated  under  the  assumption  that 
the  GB  energy  and  the  GB  mobility  are  isotropic.  The  orientation 
angle  6  assigned  to  each  grain  is  thus  expressed  by 


Fig.  2  shows  the  microstructural  features  considered  in  this 
model.  Eq.  (4)  can  be  rewritten  in  a  formulation  using  the  aniso¬ 
tropic  diffusion  tensor  in  order  to  represent  Li  diffusion  within 
a  LiCo02  grain: 


9c 

dt 
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D  = 


Dxx  Dxy 
DyX  Dyy 


DseifSin  6  +  f}self_caxiscos2^  (j^self  ^selLcaxis) 

(^self  —  ^self.caxis) Sin# COS#  ^selfcos2^  +  ^self_caxissin  ^ 

(9) 


where  Dseif_c  axis  is  the  self-diffusion  coefficient  along  the  crystal¬ 
lographic  c-axis  direction  and  6  is  the  angle  between  the  global  Li 
transport  direction  and  the  crystallographic  basal  plane  of  the 
hexagonal  LiCoCb.  The  GB  is  assumed  to  be  a  thin  layer  between 
neighboring  grains  [28].  In  the  2D  GB  model,  a  change  in  the  Li 
concentration  in  the  GB  is  calculated  by  the  following  equations, 


Layered 
Structured 
Crystallite 


Li  Transport 
Direction 


Grain 

Boundary 

Layer 


Fig.  2.  Schematic  illustration  of  microstructural  characteristics  treated  in  this 
simulation. 


6  =  mM,  (11) 

where  M  is  defined  as  180°  divided  by  the  total  number  of  grains 
(Ngb).  The  variable  m,  which  is  an  integer  between  zero  and  Ngb,  is 
randomly  assigned  to  each  crystal  grain.  When  6  is  zero  or  180°,  the 
c-axis  of  the  hexagonal  LiCoCb  crystallite  is  perpendicular  to  the 
global  Li  transport  direction.  The  diameter  of  an  irregularly  shaped 
grain  was  calculated  as  that  of  the  circle  with  the  equivalent  area.  In 
the  following  section,  the  grain  size  in  the  polycrystalline  model  is 
represented  by  the  ratio  of  the  mean  grain  diameter  to  the  length  of 
the  simulation  region  (L). 


2.4.  Numerical  calculation 


2D  numerical  calculations  were  performed  by  using  the  finite 
volume  method  with  the  explicit  Euler  method.  The  boundary 
conditions  are  listed  in  Table  1.  The  temperature  was  300.15  K,  the 
grid  spacing  was  0.02  pm,  and  the  simulation  region  was  a  square 
with  a  size  of  4.0  x  4.0  pm2.  Specifically,  a  finer  grid  spacing  near 
the  GBs  was  used  so  that  the  GB  surface  could  be  well  represented. 
The  parameters  used  in  the  numerical  calculation  are  listed  in 
Table  2.  Taking  into  account  the  experimental  Li  conductivities 
obtained  using  the  PITT  method  for  thin  films  with  (003)  and  (104) 
preferred  orientations  [10],  Dseif_c  axis  is  assumed  to  be  100  times 
smaller  than  the  Dseif  value. 

The  temporal  evolution  of  the  Li  concentration  in  the  micro¬ 
structure  induced  by  a  concentration  gradient  along  the  y-axis  was 
calculated  numerically,  as  shown  in  Fig.  4.  The  initial  concentration 
of  Li  in  the  microstructure  was  0.6,  and  the  Li  concentration  at  y  =  L 
was  fixed  to  0.61.  Once  dc/dt  is  obtained,  Dapp  of  Li  in  the  calcula¬ 
tion  region  can  be  analytically  evaluated  using  the  following 
equation  [26]: 


Dapp 


din  (dc/dt)  4L2 
dt 


if  t»L2/Dapp . 


(12) 


where  L  corresponds  to  the  macroscopic  Li  diffusion  length.  In  the 
isotropic  diffusion  system,  where  Eq.  (1 )  holds,  Dapp  corresponds  to 
the  chemical  diffusion  coefficient,  and  the  value  is  greater  than  or 
equal  to  Dseif. 
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dmeJL  =  0.05 
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Fig.  3.  Microstructures  with  randomly  oriented  grains  and  ratios  of  the  mean  grain  diameter  to  the  length  of  the  simulation  region  of  about  (a)  0.05,  (b)  0.1,  (c)  0.2,  (d)  0.3,  and  (e) 
0.4.  The  orientation  angle  of  each  grain  is  represented  by  its  gray  tone. 


When  the  electrochemical  reaction  is  explicitly  treated  with  the 
purpose  of  investigating  the  degree  of  Li  segregation  under  realistic 
discharge  conditions,  the  Butler-Volmer  equation  [13]  is  incorpo¬ 
rated  in  the  model  to  describe  the  charge-transfer  reaction 
occurring  across  the  LiCo02-electrolyte  interface  at  y  =  L.  The  total 
Li  intercalation  flux  per  unit  area  C/totai)  is  then  defined  in  the 
following  expression: 

/total  =  l  Jj(X,y  =  L) dX,  (13) 


J(x,y 


L) 


exp 


aanLiF 

RT 


(0  -  £4iCo 02) 


-  exp 


acnLiF 

RT 


(0  -  ^LiCo 02) 


(14) 


k  =  /<0COSfl(cJfx  -  cLiHcLir(cLi+r%  (15) 

where  0  is  the  electrode  potential  and  k  is  the  kinetic  rate  of  the 
electrochemical  reaction,  which  depends  on  the  angle  between  the 
basal  plane  of  LiCo02  and  the  tangent  line  of  the  interface  between 
LiCo02  and  the  electrolyte.  Here  cy  and  are  the  Li  concentrations 
in  LiCo02  and  the  electrolyte,  respectively,  and  aa  and  ac  are  the 
anodic  and  cathodic  transfer  coefficients,  respectively.  In  addition, 
ko  and  c]Tax  are  the  kinetic  rate  constant  of  the  electrochemical 
reaction  and  the  maximum  Li  concentration  in  LiCo02.  The  values 
of  these  parameters  used  for  the  simulations  are  also  listed  in 
Table  2.  The  equilibrium  potential  function  of  LiCo02  is  the  same  as 
in  Eq.  (6).  The  current  density  (/)  for  this  reaction  is  expressed  as 
I  =  7totai/f\  and  the  value  is  fixed  at  0.276  mA  cm-2.  To  obtain  the 
converged  0  value  for  each  individual  time  step,  the  Newton- 
Raphson  method  is  employed.  In  addition,  the  potential  distribu¬ 
tion  in  LiCo02  and  the  electrolyte  and  the  mass  transport  of  Li  ions 


Table  1 

Boundary  conditions. 


x  =  0  and  x  =  l 
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II 
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<< 

II 

t-i 

If 

II 
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II 

n-Vc  =  0 

c  =  constant  or  Jtotai  =  constant 

in  the  electrolyte  are  calculated  by  using  the  conventional  formu¬ 
lation  given  in  Ref.  [13]. 


3.  Results  and  discussion 


The  dependence  of  Dapp  on  the  grain  size  and  on  the  manner  in 
which  the  crystal  orientation  was  assigned  to  each  grain  was 
investigated.  We  defined  the  mean  orientation  angle  (0me an)  by  the 
following  equation: 


ny  —  ny  max  71*  ==UX  max 

$mean  —  9°-  E  E  I90  -6(nx,n y)\/(nx  _max  x  tly  max)  ? 

ny  =  1  nx  =  1 

(16) 


where  nx  max  and  nymax  mean  the  numbers  of  grid  points  along  the 
x-axis  and  y-axis,  respectively,  in  the  numerical  simulation  region. 
For  this  assignment  method,  (Wan  approaches  45°  as  the  number  of 
grains  increases.  The  Dapp  value  of  the  single-crystalline  model  with 
6  =  45°  was  about  3.6  x  10-9  cm2  s_1.  Fig.  5  shows  Dapp  as  a  function 
of  the  normalized  grain  size  ( dmeanlL )  for  simulations  using  two 
different  values  of  (3.  In  the  case  of  0  =  1.0,  Li  is  assumed  to  diffuse  in 
the  GB  with  the  same  Dseif  value  as  that  in  the  inner  grain.  At 
a  glance,  all  of  the  Dapp  values  corresponding  to  the  polycrystalline 
models  are  smaller  than  those  of  the  single-crystalline  model  with 
6  =  45°.  A  decrease  in  the  grain  size  leads  to  an  increase  in  Dapp. 
When  we  compare  the  calculation  results  for  (3  =  0.01,  in  contrast 


Table  2 

Parameter  values  for  numerical  simulation. 


Parameter 

Value 

Self-diffusion  coefficient  of  Li,  Dseif  (in  the 

1  x  1CT9  [8] 

direction  along  the  basal  plane)  (cm2  s  1) 

£>seif_c  axis  (in  c-axis  direction)  (cm2  s-1) 

1  x  lCT11 

Grain  boundary  width,  h  (nm) 

10.0 

Gradient  energy  coefficient,  k  (J  m2  moE1) 

-4.0  x  10"14 

Kinetic  rate  constant  of  the  electrochemical 

2.6  x  10“6 

reaction,  k0  (cm5/2  mol  1/2  s  1) 

Anodic  transfer  coefficient,  aa 

0.5 

Cathodic  transfer  coefficient,  ac 

0.5 

Maximum  Li  concentration,  cjjax  (mol  cm-3) 

0.051555 
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Time  (s) 

Fig.  4.  Temporal  evolution  of  the  mean  Li  concentration  in  the  simulation  region.  The 
concentration  gradually  becomes  close  to  the  concentration  supplied  at  the  position  of 
y  =  L.  The  value  of  d  In  (dc/dt)/dt  becomes  constant  at  times  larger  than  L2/Dapp. 


with  those  for  0  =  1.0,  a  decrease  in  the  grain  size  did  not  directly 
lead  to  an  increase  in  Dapp.  In  this  case,  the  GB  seems  to  have 
a  blocking  character  for  Li  diffusion.  The  Li  transport  involving  the 
GB  network  is  notably  affected  for  a  large  value  of  /3;  the  inverse 
proportionality  between  Dapp  and  the  grain  size  is  mainly  due  to  the 
Li  diffusion  process  utilizing  a  kind  of  a  channeling  through  the  GB 
network. 

Another  important  feature  shown  in  Fig.  5  is  that  there  are 
remarkable  variations  in  Dapp  even  when  dmean  is  the  same.  These 
variations  uniquely  depend  on  the  assignment  of  the  orientation 
angle  to  each  grain.  Dapp  is  plotted  against  dmean  in  Fig.  6.  The 
symbols  corresponding  to  dmeanIL  values  of  0.3  or  0.4  show  that 
Dapp  slightly  decreases  as  0mean  increases.  Since  a  6  value  of  90° 
corresponds  to  a  crystallographic  c-axis  parallel  to  the  Li  transport 
direction,  this  tendency  is  reasonable.  Flowever,  since  the  variation 
of  Dapp  is  still  wide  even  with  the  same  0mean  value,  it  is  hard  to 
estimate  the  value  of  Dapp  from  only  the  value  of  0mean.  In  fact,  the 
influence  of  the  connectivity  of  the  conduction  path  between  the 
crystal  grains  is  not  represented  by  the  6mean  value.  To  address  this, 
the  area-averaged  value  was  evaluated  using  the  following 


30 


40 


50 

0mean(°  ) 


60 


70 


Fig.  6.  Relation  between  the  apparent  diffusion  coefficient  and  the  mean  value  of  the 
orientation  angle  0mean  estimated  using  Eq.  (16)  for  (a)  /?  =  0.01  and  (b)  /?  =  1.0. 


equation  to  take  into  account  the  relative  orientation  between 
neighboring  grains: 


Fig.  5.  Relation  between  the  apparent  diffusion  coefficient  and  the  grain  size.  The  Dapp 
distribution  for  each  grain  size  and  (3  value  originates  from  differing  ways  of  assigning 
the  crystal  orientation  for  each  grain. 


fry  —  fry  .max 

^relative  *  90  — 

fry  =  1 


J2  {9O-0(nx,ny)} 


nx=l 


/  (hx.n 


:xn 


(y.max  )  • 


(17) 


In  this  case,  the  summation  takes  into  account  the  different  signs 
of  6  instead  of  the  absolute  value.  In  Fig.  7,  Dapp  is  plotted  against  the 
mean  value  of  the  relative  orientation  angle  (^relative)-  As  shown  in 
Fig.  7,  ^relative  shows  a  clear  correlation  with  Dapp.  Large  ^relative  values 
mean  that  most  interfaces  parallel  to  they-axis  have  such  a  character 
that  the  signs  of  the  angles  for  the  two  neighboring  grains  that 
sandwich  the  interface  are  opposite.  When  all  the  GBs  parallel  to  the 
y-axis  are  so-called  twin  boundaries,  the  Relative  value  reaches  90°.  In 
this  case,  the  flux  of  Li  is  restricted  by  the  region  near  the  GBs 
because  of  the  disconnectedness  of  the  2D  accessible  path  for  Li 
diffusion.  The  good  correlation  between  ^relative  and  Dapp  implies 
that  the  characteristics  of  the  GBs  parallel  to  the  Li  transport  direc¬ 
tion  significantly  affect  the  Li  diffusivity.  Recently,  experimental 
observation  has  revealed  that  the  high-coincidence  twin  boundary 
stably  exists  in  LiCo02  thin  films  [30].  Therefore,  it  is  thought  that  the 
restriction  of  Li  flux  by  the  disconnectedness  of  the  2D  accessible 


204 


S.  Yamakawa  et  al.  /  Journal  of  Power  Sources  223  (2013)  199-205 


a 


50  60  70  80  90 

^relative  (  ) 


E 

o 


2.5E-09 


2.0E-09  - 


1.5E-09  - 


1.0E-09 


^  5.0E-10 


c 

0) 

o 

*1= 

CD 

o 

O 

c 

o 

(!) 

i§ 

b 


0.0E+00 


■  •  • 


/?  =  1.0, 
^mean  f  L  — 

▲  0.05 
•  0.3 
■  0.4 


#  •*  i 


50  60  70  80 

^relative  (  ) 


90 


Fig.  7.  Relation  between  the  apparent  diffusion  coefficient  and  the  mean  value  of  the 
orientation  angle  Relative  estimated  using  Eq.  (17)  for  (a)  /?  =  0.01  and  (b)  /3  =  1.0. 


path  described  in  this  paper  occurs  in  the  actual  polycrystalline 
LiCo02  film.  Moreover,  the  Dapp  values  corresponding  to  the  dmean/L 
value  of  0.05  for  (3  =  1.0  shown  in  Fig.  7  (b)  are  somewhat  larger  than 
the  value  expected  from  the  results  for  the  dmeanlL  values  of  0.3  and 
0.4.  This  is  mainly  because  the  Li  flux  in  microstructures  with  such 
a  small  grain  size  and  such  good  conductive  GBs  is  determined  by  the 
GB  diffusion.  When  GB  diffusivity  is  poor,  e.g.,  when  /?  =  0.01,  the 
Dapp  value  corresponding  to  dmeanlL  =  0.05  is  in  good  agreement 
with  the  one  predicted  using  other  dmeanlL  values,  as  shown  in  Fig.  7 
(a).  Consequently,  both  the  GB  diffusivity  and  the  relative  orientation 
angle  between  neighboring  grains  affect  Li  diffusivity.  When  the 
grain  size  is  significantly  small,  the  Li  diffusivity  is  probably  rate- 
determined  by  the  GB  diffusion,  since  there  are  a  number  of  ways 
to  assign  the  crystal  orientation  for  each  grain  such  that  the  resultant 
^relative  values  are  almost  the  same  among  the  models  generated  for 
simulations.  On  the  other  hand,  the  diffusivity  of  a  large-grained 
structure  model  depends  considerably  on  the  relative  orientation 
angle  between  neighboring  grains,  even  when  the  Li  conductivity  of 
the  GB  is  large,  resulting  in  the  variations  of  Dapp.  When  the  electrode 
size  is  quite  large,  such  a  local  inhomogeneity  is  fully  averaged. 
However,  in  the  case  of  thin  film  electrodes,  dmean  relative  to  the 
electrode  thickness  is  not  very  large.  Therefore,  the  evaluation  of 
dmeanlL  is  useful  for  obtaining  the  appropriate  electrode  micro¬ 
structure.  Once  we  can  experimentally  evaluate  the  grain  size 


relative  to  the  electrode  size  from  a  2D  image,  we  can  roughly  esti¬ 
mate  the  possible  in-plane  variation  in  electrode  performance  in 
comparison  with  dmeanlL  squared. 

These  results  were  also  used  to  obtain  the  ratio  of  the  Dapp  of  the 
microstructure  with  randomly  oriented  grains  relative  to  that  of  the 
microstructure  in  which  the  grains  are  oriented  with  6  =  0°.  Fig.  8 
plots  this  ratio  against  the  grain  size.  The  ratio  decreased  as 
the  grain  size  increased,  and  the  ratio  varied  within  the  range  of 
0.05-0.45.  In  addition,  because  the  ratio  behaved  similarly  for  both 
=  1  and  (3  =  0.01,  it  follows  that  the  influences  of  the  GB  diffusivity 
and  ^relative  can  be  evaluated  separately.  The  experimental  value  of 
the  ionic  conductivity  obtained  using  the  electron  blocking  method 
with  a  randomly  oriented  sample  is  between  the  conductivities  for 
the  (003)  and  (104)  orientations  of  LiCoCb  [10].  Furthermore,  the 
ratio  of  the  conductivity  of  the  (104)  oriented  film  to  that  of  the 
randomly  oriented  sample  is  roughly  0.3  for  a  composition  of 
Lio.65Co02  [10].  Hence,  the  calculated  results  are  consistent  with  the 
experimental  results.  On  the  other  hand,  the  Dapp  with  =  1  is  still 
large  in  comparison  with  the  experimental  DChem  of  the  PLD  film 
measured  by  PITT,  which  varied  in  the  range  from  10~12  to 
1CT10  cm2  s-1  [9,10].  To  obtain  a  good  match  between  the  value  of 
Dapp  and  the  measured  value,  an  extremely  low  Li  GB  diffusivity  is 
suggested.  Recently,  the  results  of  first-principles  calculations  [30] 
have  shown  that  the  activation  energy  of  Li  diffusion  near  the  high- 
coincidence  twin  boundary  is  higher  than  that  of  the  crystal  inte¬ 
rior.  This  type  of  GB  seems  to  act  as  a  prevention  factor  for  Li 
diffusion.  Contrary  to  these  expectations,  recent  ESM  results  re¬ 
ported  a  high  Li  mobility  near  certain  GBs  [31].  Therefore,  we 
cannot  find  conclusive  evidence  to  clarify  the  actual  GB  diffusivity. 

Another  important  aspect  of  the  charging  or  discharging 
dynamics  is  the  Li  segregation  in  the  polycrystalline  LiCo02  thin 
film  induced  by  the  microstructural  anisotropy.  To  investigate  that, 
we  used  the  phase-field  model  in  combination  with  the  electro¬ 
chemical  model.  Fig.  9  shows  the  change  in  the  Li  concentration  in 
polycrystalline  LiCo02  during  the  discharge  process  at  a  constant 
current.  Li  tends  to  rapidly  diffuse  in  the  direction  along  the  basal 
plane.  As  shown  in  Fig.  9,  more  remarkable  Li  segregation  was 
induced  when  dmeanlL  was  larger  than  0.1.  The  simulation  results 
also  show  that  the  Li  segregation  was  enhanced  by  decreases  in  the 
GB  diffusivity.  When  the  grain  size  relative  to  the  LiCo02  thickness 
is  small,  such  as  when  dmeanIL  =  0.05,  the  existence  of  several 
possible  diffusion  paths  leads  to  a  decrease  in  the  degree  of  Li 
segregation.  These  results  suggest  that  Li  piles  up  in  one  grain  when 


Fig.  8.  Relation  between  the  grain  size  and  the  ratio  of  the  Dapp  value  of  the  micro¬ 
structure  with  randomly  oriented  grains  to  that  of  the  microstructure  in  which  the 
grains  were  oriented  with  6  =  0°.  Open  squares  and  filled  circles  denote  the  mean 
values  of  the  ratio.  The  ranges  of  variation  of  the  ratio  are  shown  as  error  bars. 
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Fig.  9.  Change  in  Li  concentration  during  a  constant-current  discharge  process.  The  Li  transport  direction  is  parallel  to  the  y-axis.  Here  cmean  indicates  the  Li  composition  in  LixCo02 
averaged  over  the  entire  simulation  region.  The  electrolyte  near  LiCo02  is  represented  by  a  blue  thin  layer  at  y  >  L. 


an  intergranular  angle  mismatch  is  large  or  when  the  GB  diffusivity 
is  poor,  while  the  assumption  of  isotropic  diffusion  is  valid  for 
strictly  limited  cases. 

4.  Conclusions 

We  have  used  the  phase-field  model  to  investigate  Li  diffusion  in 
polycrystalline  Li*Co02  randomly  oriented  on  a  nanometer  length 
scale.  The  simulation  results  demonstrated  that  Li  diffusivity  in 
anisotropic  material  was  determined  by  the  balance  between 
several  basic  characteristics  related  to  the  microstructure,  such  as 
the  GB,  the  crystallographic  orientation,  and  the  grain  size.  The 
diffusivity  of  the  small-grained  structure  is  determined  by  the  GB 
properties.  On  the  other  hand,  when  the  grain  size  is  large,  the 
variation  in  the  mean  relative  orientation  angle  between  neigh¬ 
boring  grains  increases.  Consequently,  the  diffusivity  of  the  large¬ 
grained  structure  varies  considerably  even  when  the  Li  conduc¬ 
tivity  of  the  GB  is  large.  It  is  possible  that  the  diffusivity  of 
a  randomly  oriented  and  large-grained  microstructure  will  be  one 
order  of  magnitude  smaller  than  that  of  the  6  =  0°  textured 
microstructure.  This  work  has  provided  a  clue  to  understanding  the 
Li  transport  properties  such  as  those  in  a  thin  film  electrode  for  an 
all-solid-state  battery.  Since  observable  phenomena  only  represent 
consequences  of  the  multistep  process,  it  is  difficult  to  separately 
evaluate  the  role  of  each  basic  step.  Therefore,  the  integrated 
computational  approach  that  simultaneously  simulates  the 
microstructure  and  evaluates  the  transport  properties  presented 
here  may  help  to  obtain  a  quantitative  understanding  of  each  basic 
step  of  the  Li  transport  process. 
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